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ABSTRACT 



Context. HD 141569, a triple star system, has been intensively observed and studied for its massive debris disk. Until recently, it 
was rather regarded as a gravitationally bound triple system but recent measurements of the HD 141569A radial velocity seem to 
invalidate this hypothesis. The flyby scenario has therefore to be further investigated to test its compatibility with the observed disk 
structures. 

Aims. We present a study of the flyby scenario for the HD141569 system, by considering 3 variants: a sole flyby, a flyby associated 
with one planet embedded in the disk and a flyby with two planets in the disk. We discuss the merits of each one to reproduce the 
scattered light observations of the disk. 

Methods. We first use analytical calculations to reduce the parameter space of the 2 stellar companion's orbit and then perform N- 
body numerical simulations of the flyby encounter, using symplectic integration, taking into account the gravitational influence of the 
stars and the planets on massless test particles. 

Results. The binary orbit is found to be almost fixed by the observational constraint on a edge-on plane with respect to the observers. 
If the binary has had an influence on the disk structure, it should have a passing time at the periapsis between 5000 and 8000 years 
ago and a distance at periapsis between 600 and 900 AU. It also appears that the best scenario for reproducing the disk morphology 
is a flyby with 1 planet embedded in the disk. For a 2 My planet, its orbital eccentricity must be around 0.2 while for a 8 Mj planet, it 
must be below 0.1. In the two cases, its apoapsis is about 130 AU. 

Conclusions. Although th global disk shape is reasonably well reproduced, some observed features cannot be explain by the present 
model and the likehood of the flyby event remains an issue for the scenario explored in this paper. Dynamically speaking, HD 141569 
is still a puzzling system. 
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1. Introduction 

HD 141569, a triple star system, has been intensively observed 
and studied for its circumstellar dusty disk imaged in scattered 
light in the visible and at near-infrared wavelengths. This system 
is located at 99 ± 10 pc an d its age i s estimated to 5 ± 3 Myrs 
jWeinberger et alj 2000; M erin et alj|2004) . The disk is associ- 
ated to the central star (B9.5 Ve) while the two other companions 
. (M2 and M4) form a binary. 

The resolved images of the disk, showing a complicated 
morphology, have generated a lot of discussion in the litera- 
ture about the systems dynamics (see Table [T] for a summary). 
First of all, this disk is in transition to a debris disk: although 
it is a dusty optically thin disk w ith a fractional disk lumi- 
nosity of Ldisk/L, ^ 8.4 X 10"^ dSvlvester et all [1996). the 
gas mass is not negli gible. It could in fact represent most of 
the mass according to iJonkheid et al.l (12006') (Mgas - 80 M®, 
Mdust = 2.2 Me for grains with radius sizes between 1 yum 
and 1 cm), although the spatial distribution of the gas and of 
the dust may differ. The influence of the gas on the dust dy- 
namic s has been therefore taken into account by lArdila et al] 
(l2005h . but other authors consider only classical N-Body simu- 
lations withou t gas ( Augereau &^ Pap aloizou 2004; Wyatt 2005). 
lOuillen et al.l J2005h . on the other hand, only considered the gas 
in their simulations. 
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The dynamical status of the external binary is also question- 
able: in order to reproduce the external structures of the disk, 
previous studies ha ve either considered the case of gravitational 
bound companions jAugereau & P apaloizou'2004';'O uillen et al.l 
2005) or the case of a flyby dXrdila et al. 2005). The internal 
disk structures might be explained by the gravitational pe rturba- 
tion o f an unseen planet at large distance, around 250 AU dWvattl 
|2005|) but also by alt ernative mechanisms, such as d ust migra- 
tion in a gaseous disk jTakeuchi & Artvmowicd200lb . The latter 
scenario is successful in producing annular structures, but can- 
not account for the observed non-axisymmetric of these features. 
The combined effect of a planet and external perturbers has only 
been considered i n one study, witho ut successfully reproducing 
all the structures dArdila et al.ll2005l) . 

Ne w millimeter observations of the gas disk dPent et all 
(l2005h and Augereau et al., in prep) , which give better con- 
straints on the radial velocity of the primary star, show that 
the differential velocity between the binary and the central star 
(5.8 ± 0.3 km s ') is larger than the system escape velocity (2.6 
km s"^) so that the flyby scenario seems therefore to be the 
most plausible. Therefore, we propose in this paper to exten- 
sively study the flyby scenario, as it is the least studied up to 
now. In Sect. |2] we summarise the available observations of the 
disk and the stars astrometry constraints. In order to best repro- 
duce the observations, we discuss 3 different scenarios (Fig.[TJ: 
a sole flyby (Sect.[3j, a flyby with 1 planet embedded in the disk 
(Sect, m and a flyby with 2 planets embedded in the disk (Sect. 
|5]l. In Sect.|6l we compare our approach to previous studies and 
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Table 1. Summary of recent papers on the dynamical modelling of the HD 141569system. 



Authors 



Modelled structures 



External companions Planet 



Notes 



Take uchi & Artvmowicz (2001) 
Aug ereau & Papaloizou (2004) 
Guillen et al. (2005) 
Ardila et al. (2005) 
Wvatt (2005) 



gap & two rings 

outer ring 

outer ring 
spiral arms & inner depletion 
gap & outer ring 



No 
bound 
bound 
flyby 

No 



No Dust migration due to gas friction 

No N-Body gravitational code 

No 2D hydrodynamics code 

Yes N-body code, SPH code, collisions 

Yes N-body gravitation code 





Fig. 1. Illustrations of the 3 scenarios studied in this paper: a sole flyby (left panel), a flyby + 1 planet (middle panel) and a flyby 
+ 2 planets (right panel). The stars and planet locations, are represented by large red points, and the various orbits by a thin green 
line. See the electronic edition of the Journal for a color version of this figure. 



discuss in particular the likelihood of this scenario. We finally 
summarise and conclude in Sect.|7] 



2. Available Observations 

2.1. Overview of the disk shape 

The motivation for our study is to reproduce the dust disk struc- 
tures that have been observed in scattered light. According to 
the visible, near-infra r ed and mid-infrared ob s ervations from 
lAugereau eTaTI d 19991). [Weinberger et all d 19991). iMouillet et all 



lAugereau et al. ( 19991). Weinberger et al. ( 
( l200lh . lMarsh et al.l ( l2002h and lClampin et 
ing facts appear: 



alj(l2003|)thefollow- 



- The main structure is the two rings shape of the disk. The 
two annuli peak at ~ 200 AU and ~ 325 AU from the star. 
The two bright rings show out of phase brightness asymme- 
tries of up to factors of 2.5-3 for the outer ring in the visible. 
The outer ring moreover shows a tightly-wound spiral struc- 
ture. The two annuli are separated by a darker ring or "gap", 
which is radially wide compared to the two annuli. 

- an extended diffuse emission associated with a faint spiral 
arm is present in the North-East of the disk and is detected up 
to more than 600 AU. According to the millimeter maps and 
using the hypothesis of forward scattering, this spiral would 
be a trailing structure with respect to the disk rotation. An 
other spiral arm is possibly observed, pointing toward the 
companions. 

- the disk brightness sharply decreases between 200 and 150 
AU from the central star until the background level is 
reached. This behaviour is suggestive of a strong depletion 
of dust inside 150 AU. 

2.2. Radial velocities & astrometry 

[Wein berger e t al.l (|200(I ) have measured the radial velocities of 
the HD 141569B and C stars, respectively - 1 .5 + 0.6 km s ' and 
-2.4 ± 0.7 km s"' by cross-correlating their spectra to that of 
standard stars. This method was not suitable for HD 141569A 



due to its lack of lines for comparison with the radial velocity 
standard. [Weinberger et al.l ([2000) also summarise the astrome- 
try of the 3 stars from 1938 to 1999^ 

Augereau et al. (in prep) and iDent et aTl ( l2005l) have mea- 
sured the double-peaked CO J=2-l spectral profile of the disk 
around HD 141569A. It is thus possible to deduce the radial ve- 
locity of this star from the CO line and to obtain a far smaller 
uncertainty on this measure. To a first approximation, assuming 
a symmetric emission in the blue-shifted and red-shifted parts of 
the disk, the radial heliocentric velocity of HD 141569A is found 
to be -7.6 + 0.3 km s"' (mean velocity in the line). Supposing 
that the star masses are 2.5, 0.5 and 0.25 Mq for A, B and C 
respectively ( I Weinberger et al.l |2000|) and that the distance be- 
tween HD141569A and the center of mass of the B-C binary is 
at least 800 AU (projected distance, [Weinberger et all ([2000')), 
the maximum escape velocity is 2.6 km s~' . As the measured 
velocity is far above this limit, the two M star companions are 
not bound to HD 141569A. In the context of dynamical pertur- 
bation of the disk by the companions, the flyby scenario seems 
therefore to be the most plausible. 



3. The flyby scenario 

3.1. Kinematic constraints 

We define a cartesian coordinate system {X, Y, Z) where X points 
toward the north, Y toward the east and Z toward the Earth and 
where the origin is HD141569A. In this referential frame, the po- 
sition of the center of mass of the B-C binary will be described 
by the coordinates (Xf,, Yi,,Zi,) and its velocity by V/, , V/, , Vh.. 
We already know the present Xi,, Yt (right ascension and dec- 
lination) and Vh^ (radial velocity). But the hyperbolic orbit of 
the binary around the primary is also defined by the 6 orbital 
elements (q, e, i, o), D., u) where q is the pericenter distance, e 
the eccentricity, / the inclination with respect to the sky plane 
iOXY), (jj the periapsis argument, O the longitude of the ascend- 
ing node and u the eccentric anomaly at the observing time. The 
transformation from cartesian coordinates to orbitals elements, 
or the inverse transformation, is defined by a system of 6 non- 
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linear equations from classical keplerian formalism. It is thus 
possible to look for orbital element sets verifying the 3 observ- 
ables (Xh, Yb, Vzb)- But with 6 unknown and only 3 constraints 
the equation system is underconstrained. One has therefore to 
consider 3 orbital elements as free parameters in order to solve 
the system for the 3 other orbital elements. This system is de- 
fined by the following equations, for an hyperbolic orbit: 

Xf, cos Q. - Yh sinQ. = a cos / ^ Ve^ - 1 cos lo sinh u 

+ (e - cosh u) sin cS) (1) 

Xb sin Q + Ifo cos il = a Ve^ - 1 sin lo sinh u 

+ (e - cosh u) cos oj) (2) 
Vfc. Vfl cosh M - 1) = VGM sin / (sinh u sin oj 

- Ve^ - 1 cosh M cos (3) 

where G is the gravitational constant and M the sum of the 3 
star masses. For this study, we choose to consider (/, Q., u) as free 
parameters and {q,e,io) as unknowns. Once the 3 parameters 
{i, n, m) are fixed to an arbitrary set of values, the non-linear sys- 
tem generally admits none or one solution for the triplet {q, e, a>). 
For our problem, many triplets of the parameters (/, Q, u) give a 
solution: for each (/, M)pair , at least one value of Q. gives a valid 
solution for (q, e, oj). The results are thus not significant because 
all the scenarios can be compatible with the observations: orbital 
plane perpendicular to the line of sight, or instead very inclined 
orbits, more or less eccentric hyperbolic orbits ... 

However, an other constraint exists, not yet used in our 
method: th e observed proper motio n of the stars between 1938 
and 1998 dWeinberger et al.l r200(]h . and the associate velocity 
Vskv = 1.1 ± 0.7 km s"^ The uncertainty on this observable is 
quite large and we consider it as an upper limit: all the solutions 
that predict a proper motion less than twice the observed proper 
motion are valid: 

+ Vl < 2V,ky (4) 

This is enough to reduce significantly the number of solu- 
tions, as showed in Fig. |2l On this plot, the two axes are the 
eccentric anomaly u and the inclination /, while the colors cor- 
respond to the density of valid solutions (i.e. fulfilling the re- 
lation 1 to 4) for the longitude of the ascending node Q. in the 
range [0, 27r]. The parameter space region where the valid so- 
lutions exist is indeed very much reduced by using the proper 
motion constraint as we find, with an inclination / between 75° 
and 105° and an eccentric anomaly u between -5 and 5 radians. 
The fourth constraint has therefore eliminated all the solutions 
with lower inclination and larger proper motion. 

The disk inclin ation is 55° ± 1.2° from pole-on 
(iMouillet et al.l l200lh . or rather 180° - 55° in order to 
have the disk rotation jGoto et al.l 12006) consistent with our 
convention (i.e. ; < 90° is an anticlockwise orbit on the sky). 
The binary orbits with an inclination below (resp. above) 90° 
are thus retrograde (resp. prograde) with respect of the disk 
rotation. The relative inclination between the orbital plane of 
the binary and that of the disk can be calculated by the following 
equation: 

cos ii-ei - COS i COS i^iisii + sin i sin ; j,-.5t cos (Q - Q.disk) (5) 

with Q.disk, the position angle of the disk (357° + 2°) and /, Q 
orbital elements of the binary orbit. 
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Fig. 2. Density of valid orbital configurations in the parameter 
space (;, u, Q.). The two axis are the eccentric anomaly u and the 
inclination /, while the colors correspond to the density of the 
valid solutions for the longitude of the ascending node D. in the 
range [0, 27r]. The small white box shows the final valid param- 
eter space when the disk geometry constraints are also taken in 
account (Sect. |X2] i. See the electronic edition of the Journal for a 
color version of this figure. 



3.2. Disk geometry constraints 

At this point, only the kinematic constraints have been used to 
calculate the valid orbital configurations for the binary. However, 
the encounter geometry is also constrained by the observed 
structures in the disk if one assumes that the current disk shape 
is a consequence of the gravitational perturbations of the disk 
by the binary. The comparison study of numerical simulations 
of the flyby scenario to the observations therefore reduces even 
more the range of valid orbital configurations as detailed below. 

3.2.1. Numerical model 

We consider in our simulations a hierarchical system consisting 
of a central star associated with a debris disk and an external bi- 
nary. The trajectory of this binary can be a bound orbit around 
the primary star or a flyby (hyperbolic orbit). We address the 
case of a disk consisting of large particles (planetesimals), which 
are insensitive to pressure forces (radiation, stellar wind or gas 
pressure). Importantly, we also do not take into account the grav- 
itational interactions between planetesimals as they are negligi- 
ble, nor mutual collisions. Dynamically speaking, the planetes- 
imals are thus considered as test particles. In Sect. |4] and |5] one 
or several planets will be added to the simulations, orbiting the 
central star in the orbital plane of the disk. The central star is 
assumed to have a mass of 2.5 Mq and the low mass compan- 
ions companions 0.5 and 0.25 (Weinberger et al. 2000i). 
The initial disk consists of 100000 test particles with a surface 
density distribution proportional to r ' . 

To perform our simulations we u se the symplectic package 
HJS (Beust 2003), a SWIFT variant ( Wisdom & Holman|[T99ll: 
iLevison & Duncan, 199 4) for hierarchical systems. It allows us 
to integrate accurately and fastly the motion of the disk parti- 
cles, the relative motion of the two companions and their orbit 
around the central star, although they are on different timescales. 
Numerous scenarios can be therefore tested and studied in a rea- 
sonable CPU time. 
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3.2.2. Results 

Among the solutions consistent with the kinematic constraints 
(Sect 13. It . the periapsis of the orbit, do, (i.e. the closest ap- 
proach) can range between 10 and several thousands of AU. If 
the flyby is however too close to the star, the disk can be totally 
dissipated or, at least, can show a spiral structure which strongly 
differs from the observations (Fig.[3]top left panel). We numer- 
ically explored a broad range of do values consistent with the 
kinematic constraints of Sect. 13.11 According to the simulations, 
we find that the distance at periapsis, do must be above about 
600 AU. An other useful constraint is the time span between the 
periapsis passage and the present time, T = f - fo- If ^ is too 
small compared to the keplerian period at the outer disk edge, 
the spiral structure does not have the time to develop (Fig.[3]top 
right panel). The simulations give an estimate of at least 5 000 
years before present for the periapsis passage. All the well de- 
fined structures in the disk develop after the periapsis passage, 
never before, which eliminates all the configurations with m < 0, 
because the binary is still too far to generate any perturbations in 
the disk (Fig.[3]bottom left panel). 

Few orbital configurations finally verify the system given by 
the equations[T]to|4]plus the additional constraints do > 600 AU 
and T = t - to > 5 000 yrs. Valid solutions are only possible 
for / between 85° and 95° and u x 3 and share the following 
properties: 

5000 yrs <t- to < 8000 yrs 
600AU<do <900AU 

e a: 10 

Thanks to Eq.|5] it is also possible to determinate that: 

35° < < 145° 

Considering these solutions, the binary has at present time 
already reached its velocity at infinity. Within this range of solu- 
tions, the initial size of the disk can also be better constrain by 
our simulations in the range between 400 - 450 AU. If it is larger 
(Fig. |3]bottom right panel), the flyby produces large spiral arms 
which are too extended while if the disk size is smaller, the disk 
does not fit the observed size. We do not address here the ques- 
tion of the origin of the disk size, nor the shap e of the external 
disk edge. But collisional evolution of the disk (iThebault & Wul 
12008 ) or gravitationnal truncation by the binary in a past bound 
state of the system (see Sect. 16. U can both be used to explain 
them. 

Overall, this scenario gives satisfactory results for the global 
shape of the disk and reduces significantly the range of possible 
orbital configurations for the companions. This is illustrated in 
Fig. |4] which shows a disk of correct size with a outer trailing 
spiral. The inner disk structure is nevertheless not properly re- 
produced and the next step is thus to reproduce more accurately 
the two ring-li ke structure at distances of about 200 and 3 25 AU 
from the star jMouillet et al.ll2001l: IClampin et aLll2003l) . It ap- 
pears however that none of the valid solutions for the flyby is 
able to well generate the two rings. In the case of retrograde 
orbit (85° < / < 90°), the flyby generates a clear two arms spi- 
ral structure (Fig. ID left panel). In the ca se of p rograde orbit 
(90° < / < 95°), as studied for instance by|PMzner| (|2003|) . the 
two arms are more wound and more extended but it is still too 
different from the observations (Fig. |4] right panel). In the two 
cases, the flyby does not perturb the inner part of the disk, below 
200 AU. It thus cannot explain the inner ring, nor the observed 
depletion around 150 AU. It means that we need to add some 
other hypothesis to the model. 




Fig. 4. Best solutions for the flyby scenario on a retrograde orbit 
(left panel) and on a pr ograde (right panel orbit, compared to 
the observations from Moui llet et al.l (12001 ) (middle panel). A 
scattering anisotropic factor g of 0.2 is assumed to mimic real 
scattered light observations. 

4. The flyby + planet scenario 

We can conclude from Sect. [3] that the sole presence of the 
two star companions is not enough to well explain the observed 
disk properties summarised in Sect.|2] especially the double-ring 
structure. The flyby can, at least, constrain the initial disk size 
between 400 and 450 AU but it is necessary to use other hy- 
pothesis to improve the model and to try to better fit the disk 
geometry. As the observations show a strong depletion in the 
dust surface density below 150 AU, a plausible hypothesis is the 
presence of an unseen planet which generates this gap. 

In this section, we consider the case where we add a planet 
within 150 AU in our numerical models. The planet is supposed 
to be coplanar with the disk and its main parameters are its mass, 
its orbital semi-major axis and its orbital eccentricity. The sim- 
ulations take into account the mutual gravitational interactions 
between the planet and the 3 stars, as well as the gravitational 
perturbation of the disk by the planet but not the feedback of the 
disk on the planet: the planetesimals are still considered as test 
particles. The system consisting of the central star, the disk and 
the planet is first evolved during 5 Myrs (the assumed age of the 
system), then the two star companions are introduced to simulate 
the flyby with one of the valid set of orbital parameters found in 
Sect. [31 

4. 1 . The inner ring 

To constrain the planet parameters, we co mpare the simulation s 
to the disk optical thickness derived from iMarsh et al.l ( l2002h . 
The authors h ave combined their ow n observations in the mid- 
infrared and Weinberg er et al.l (1 19991) observations in the near- 
infrared (Fig. Is) to estimate the radial profiles of the vertical dust 
optical depth at 1.1 fim. The observations in the mid-infrared 
cover the inner region of the disk, up to 130 AU while the ob- 
servations in the near-infrared cover the region between 130 and 
500 AU. The radial profile of the optical depth at 1 . 1 /im is de- 
duced, for the mid-infrared data, from the optical depth at 12.5 
fim and assuming the opacity scales inversely with wavelength. 

The profiles show a ring which peaks around 200 AU, a de- 
pletion at 250 AU and an even more depleted region below 150 
AU. We use in the following these data to constrain the mass and 
the orbital elements of the planet, by trying to reproduce some 
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Fig. 3. Examples of parameters for the flyby scenario inconsistent with the observations: the periapsis is too close to the disk, at 170 
AU, (top left panel) the crossing time at periapsis is too recent, 700 years, (top right), the binary is still not yet at periapsis (bottom 
left) and the disk is too large (bottom right). These figures plot the intensity, in logarithm scale, of scattered light by the planetesimal 
spatial distribution obtained with our simulations. The scattering anisotropic factor is 0.2. The time is relative to the crossing time 
at periapsis, fo. 



key characteristics of the azimuthal averaged radial density pro- 
file of the dust: 

- a peak between 190 and 210 AU 

- a ratio t peak I tiso of about 1 .3 between the peak surface den- 
sity at ~ 200 AU and the depletion at ~ 250 AU 

- a ratio Tpeak/Ti^o of about 2.7 between the peak at ~ 200 AU 
and the inner depleted region at 150 AU. The uncertainty 
on this ratio is nevertheless rather large because of the large 
uncertainties in the optical depth between 100 and 150 AU. 

The "1 planet + flyby" scenario is tested for 3 different planet 
masses: 0.2, 2 and 8 My. For each planet mass, 3 apoapsis are 
considered (130, 150 and 180 AU) with 4 diff'erent values of the 
eccentricity: 0, 0.05, 0.1 and 0.2. This gives a total of 36 simula- 
tions and their results are summarised in Tables |2] E] and |4] The 
peak position (Table |2]i indicates the position of the largest over- 
density if several are present. Tpeak /tiso is defined as the ratio 
between the dust density at the peak position and the maximum 
density between 100 and 150 AU. t peak 1^150 is defined as the ra- 
tio between the dust density at the peak position and the density 
around 250 AU. 

According to these simulations, a 0.2 My planet cannot re- 
produce the observational data: it does not clear out enough the 
inner part of the disk. For higher mass planets, the best solu- 
tion is obtained for an apoapsis of 130 AU. For larger values of 
apoapsis, the peak position is shifted toward too large distance 
from the star, or the planet can even prevent the formation of an 



over-density. For a 2 My planet, the best fits are obtained for ec- 
centricities around 0.2, while for a 8 My planet, it must be below 
0.1. The best solution between these two planets is for a 2 My 
planet. 

With an initial surface density distribution proportional to 
r"', Fig.|5]shows that there is a lack of particles in the outer part 
of the disk compared to the optical depth profile al A - 1.1 yum. 
We therefore adjust the initial surface density to obtain the same 
maximum at the outer peak location, near 330 AU. We find that 
the best fit is obtained for an initial surface density distribution 
proportional to r^^-^ and that the outer peak position in our simu- 
lation is consistent with the observations. However, there is still a 
lack of particles between the two rings and outside the outer ring. 
This might reflect a limitation of our model as we only model 
the dynamics of planetesimals. One can expect than smaller dust 
particles, sensitive to radiation pressure and generated by col- 
lisions, will migrate outside the rings and flatten the profiles, 
hence better reproducing the optical depth profile. 

4.2. The global structure 

The ratio between the two peaks in the optical depth profile de- 
pends of the observing wavelength, as the optical depth at the 
outer peak position is larger than the optical depth at the inner 
peak position at 0.5 fim ( Ardila et al. 2005), while it is the oppo- 
site of Fig.|5](/i =1.1 fim). Our constraint on the initial surface 
density distribution should therefore be taken with care, but our 
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Table 2. Summary of the results for the peak position (in AU) in the planetesimal density profile, for all the simulations, as a 
function of the planet mass, eccentricity and apoapsis. "-" indicates the simulations where no peak is observable or where the peak 
is below 170 AU. t indicates the simulations where multiple peaks appear. In this case, the largest peak is used to calculate the peak 
position, T peak 1^X50 and TpeaklTtso- The numbers boldface correspond to the best solutions. 



Planet mass 
Planet ecc. 



apoapsis of 180 AU 
apoapsis of 150 AU 
apoapsis of 130 AU 
apoapsis of 120 AU 



0.2 My 



0.05 



0.1 



0.2 



2 Ml 



0.05 



0.1 



0.2 



8Mj 



0.05 0.1 0.2 



210 
220t 



210 



210 



230 
220 



230 
2I0t 
220t 



230 



200 



240 
210 
190 



220 220 230 220 
220 220 230 230 



Table 3. Same as Table|2l but for the ratio TpcaklTuo between the density at peak position and the maximum density between 100 
and 150 AU. 



Planet mass 
Planet ecc. 



0.2 M, 



0.05 



0.1 



0.2 



2M, 



0.05 



0.1 



0.2 



8 My 



0.05 



0.1 



0.2 



apoapsis of 180 AU 
apoapsis of 150 AU 
apoapsis of 130 AU 
apoapsis of 120 AU 



0.5 
0.5t 



0.5 



0.5 



0.8 
1.1 



0.6t 
2.4t 



0.4 
1.6 

1.3 



3.0 
2.3 
21 



3.0 
2.7 



5.0 
4.7 



6.0 
4.2 



8.0 
4.5 



Table 4. Same as in Table|2] but for the ratio Tpeak/T250 between the density at peak position and the density in the gap around 250 
AU. 



Planet mass 
Planet ecc. 



0.2 M, 



0.05 0.1 0.2 



2M, 



0.05 0.1 0.2 



8 My 



0.05 0.1 0.2 



apoapsis of 180 AU 
apoapsis of 150 AU 
apoapsis of 130 AU 
apoapsis of 120 AU 



1.5 
1.7t 



1.7 



1.3 



1.8 
2.0 



1.1 

1.4t 1.5 2.0 

1.7t - 2.0 

2.2 2.1 



1.9 
1.8 



1.9 
1.7 



1.7 
1.8 



1.3 
1.4 



Mid-infrared data i — i — 
Near-infrared data '-— x— 
Simulation ^ 

1 Simulation witli constrained surface density □ 



to 



^ □ 



200 300 
r[AUl 



Fig. 5. Estimated profiles of dust optical depth ( normal to the 
disk plane) by combini ng mid-infrared d ata ( Marsh et ani2002h 
and near-infrared data (IWeinberger et al.llT999) . compared to the 
simulation output of the best scenario, a 2 My planet with an 
eccentricity of 0.2. 



model is nevertheless very successful to predict the position of 
the outer peak in the optical depth profile. The planet parameters 
are strongly constrained by the inner ring, but our simulations 
show that the planet also perturb the outer parts of the disk as 
shown in Figs.|5]and|6] 



In less than 5 Myrs, the planet has time enough to develop a 
strong, one-armed, anticlockwise wound spiral, due to the differ- 
ential precession rate of the planetesimal orbits as explained by 
lAugereau & Papalo izou (2004). The flyby does thus not occur 
with an initially smooth disk and this has strong consequences 
on the final shape of the structures produced by the flyby. For 
instance, in the case where the binary is on a retrograde orbit, 
we previously obtained a regular two armed spiral. But now, the 
planet not only produces the inner ring, but can also make disap- 
pear, or at least perturb, a significant fraction of one of the two 
spiral arms. As shown in Fig. |6] the disk morphology between 
200 and 300 AU slightly depends of the direction of the planet 
periapsis with respect to the direction of the binary periapsis, 
but the differences remain too small to constrain the orientation 
of the planet orbit. The results are thus not perfect because the 
simulations do not still reproduce exactly the outer ring but they 
clearly show that we need a complete description of the system 
to properly account for it. It is not satisfactory to consider solely 
the external companions or the inner planet. Combining both, we 
clearly achieve a much better fit to the observations. 



5. The flyby + 2 planets scenario 

As the previous scenario is still not fully satisfactory because 
it still does not generate two rings clearly separated by a large 
gap, we have also investigated an alternative with two planets 
in the disk. Wyatt (2005) has already considered the case of a 
single planet at about 250 AU to explain the outer ring but with- 
out considering the perturbation due to the companion stars. We 
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Fig. 6. Same as in Fig.|4l but for the best solution for the flyby + planet scenario in the case of a retrograde orbit of the binary (third 
panel). The red arrow indicates the direction of the planet periapsis. The planet mass is equal to 2 My, it apoapsis is 130 AU and 
its eccentricity 0.2. The first panel show the disk before the flyby but after 5 Myrs of gravitational perturbation by the planet. The 
fourth panel shows an alternative configuration where the initial direction of the planet periapsis is turned by 180°. 

be considered: a 0.2 Mj planet evolving over 5 Myrs on an orbit 
with a semi-major axis of 250 AU and an eccentricity of 0.05, 
and a 2 Mj planet perturbing the disk over 0.5 Myrs and on the 
same orbit. We add thus another planet in the internal part of the 
disk, following the best solutions found in the previous section, 
namely a 2 My planet on an orbit with a semi-major axis of 1 10 
AU and an eccentricity of 0.2 or a 8 My planet on an orbit with a 
semi-major axis of 1 15 AU and an eccentricity of 0.05 (Fig.|7]i. 

Like in Sect. |4] the estimated radial profiles of the dust op- 
tical depth obtained with these simulations are compared to the 
observed optical depth at 1.1 yum (Fig.lSj. Two major informa- 
tions are therefore given by this method. First, and independently 
of the mass of the inner planet, a 2 My at 250 AU generates a 
too deep gap in the disk compared to the observations (scenario 
B and C of Fig. [SJ. Contrary to scenario "flyby + 1 planet", 
one cannot expect that the optical depth profile would be flat- 
tened by smaller dust particles sensitive to radiation pressure be- 
cause the particles will be quickly removed from the ga p, due 
to clo se encounter with this planet orbiting at 250 AU. IWvattI 
(l2005b constrains the planet parameters considering the location 
and the winding of the spiral but there is a degeneracy between 
the planet mass and the time the planet is orbiting around the 
star. If one adds the constrain given by the observations on the 
depth of the gap around 250 AU, only the scenarios with smaller 
planet masses (0.2 My) orbiting for a long time (5 Myrs) are sat- 
isfactory (scenario A and D of Fig.|8]l). Second, the simulations 
show that it is difficult to obtain an inner ring massive enough 
between the two planets, around 200 AU. Due to the perturba- 
tions of these two objects, a lot of planetesimals are ejected from 
the disk. We have therefore tried other simulations with a smaller 
mass (0.2 My, scenario D) for the internal planet in order to less 
perturb the disk. However, the results are still not satisfactory 
because the peak position is not at the correct location (150 AU 
instead of 200 AU) and because, with a low mass planet, there 
are too much planetesimals below 100 AU. 

A planet located at 250 AU from the primary star is a in- 
teresting alternative to the companions in order to explain the 








1 





Fig. 7. Same as in Fig. |4], but for the flyby + 2 planets scenario. 
The inner planet mass is equal to 2 My, it apoapsis is 130 AU and 
its eccentricity 0.2. The outer planet mass is 0.2 My, it apoapsis 
is 250 AU and its eccentricity 0.05. 



therefore first investigate the consequence of both the compan- 
ions and the planet on the disk geometry. Using the solutions 
found in Sect. |3] it appears that adding the gravitational influ- 
ence of the two companion stars does not improve the structure 
generated by the external planet but, on contrary, disrupts them 
because the planetesimals orbits are too destabilised by the per- 
turbations of a planet and that of the two companions. For this 
scenario, the constraints used in Sect. |3.2| ai'e no longer valid and 
the correct orbital elements for the binary orbit are those which 
prevent interactions between the disk and the two M stars: they 
do not still reach the periapsis location (m < 0) or the periapsis 
itself is located at more than 1500 AU from the primary star 

Such a scenario cannot, however, reproduce the inner deple- 
tion below 150 AU; another pl anet is thus n ecessary like in the 
previous section. According to lWvattI (l2005l) . two scenarios can 
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Fig. 8. Optical depth obtained by simulations with 4 different 
scenarios for the mass of the 2 planets, to compare to Fig|5] 



outer ring without the gravitational perturbation of the 2 stel- 
lar companions, but our work shows that combining it with an 
other planet around 120 AU to explain the inner ring and the 
dust depletion within 150 AU does not provide satisfactory re- 
sults. Even if it is still not perfect, the scenario "flyby + 1 planet" 
is the one which gives the best results among the three studied 
scenarios in this paper. 

6. Discussion 

6. 1 . Probability of tiie encounter 

The main drawback to the flyby scenario is that, according to 
our models, the binary passed to the periapsis less than 8000 
years ago. The probability of observing such an event is low, 
making the flyby scenario statistically unlikely. On the other 
hand, the similar physical properties between the central star 
and the low mass companions suggest a common origin. The 
3 stars have been independently estim ated to be around 5 Myrs 
([Weinberger et all 1200(31: iMerfn et al.|[2004 ) and the diff'erential 
radial velocity between the HD 141569A and the binary is only 
a fraction of the total speed of the primary star. But the observed 
radial velocities clearly show that the binary is on an escape or- 
bit and the open spiral structure observed in the disk is also in 
favour of a flyby scenario, as trailing structures are created by 
tidal effects (Augereau et al., in prep.). 

One can however reconciliate these positions with an hybrid 
scenario. HD 141569 could have been, in the past, a triple bound 
system until an other star of its vicinity approached and desta- 
bilised the system, putting the binary on its current orbit. In this 
case, we would be observing the last periapsis passage of the 
binary and the consequences of a 4 star encounter This sce- 
nario could explain the observed velocities, the common phys- 
ical properties and the disk size, and it is more plausible that a 
sole flyby of the binary. The gravitationally bound scenario stud- 
ied by^ugereau & Papaloizou (2004) and Ouillen et al. (2005) 
can possibly be used to constraint the orbit of the binary before 
the encounter: typically an orbital period of several hundreds of 
thousands years and a apoapsis of almost 10000 AU. It appears 
more probable that a fourth star destabilises such an orbit that the 
binary randomly cross the p ath of the primary star at less than 
1000 AU. The study done bv lScallv & Clarkd ( l200lh show that, 
in the Orion Nebula Cluster, almost all the stars undergo at least 



one encounter betw een 1000 and 10000 AU. On the other hand, 
Aarnio et all (|2008[) have found that HD 141569 formed in rela- 
tive isolation, tens of parsecs away from the recent sites of star 
formation in the Ophiucus-Scorpius-Centaurus region. The like- 
hood of such a scenario has thus to be better estimated. Its dy- 
namical modelling has also to be investigated by this is beyond 
the scope of this paper. Theoretical works are however already 
abundant in the literature (see, for example, the s eries discussing 
many as pects of the three-bo dy scattering from iHut & BahcaUl 
([T983h to lHeggie etal! (199^ ). 

6.2. Comparison witli previous works 

Other theoretical studies of the HD 141569system can be divided 
in 3 groups with respect to the reproduced disk structures: the 
outer ring, the gap between the two rings or the extended spiral 
arrnS; 

lAugereau & Papaloizoul (l2004 and lOuiUen et all (l2005h 
have both studied the outer rin ^ of the disk, but using d ifferent 
tools: a pure N-body code for Augereau & Papaloizou (20041) 
and a 2D hydrodynamics code for Ouillen et al. (2005) assum- 
ing in the latter case that the gas and the dust are fully mixed. 
Their main issue compared to our work is that they have con- 
sidered bound orbits for the binary which is incompatible with 
the new measure ment of radial veloc ities for the primary star 
Interestingly, Augereau & Papaloizoul 4200 4) well reproduce the 
outer ring but not the open spiral arms structure while it is the 
contrary for our work. One can therefore ask if the scenario with 
4 stars can reconciliate these two approaches. This has to be in- 
vestigated in a fu ture paper. 

IWvattI ( l2005h and lTakeuchi & Artvmowicd ( 1200 lb have con- 
sidered the opening of the gap between the two rings wit h very 
different mechanisms: a planet embedded in the disk dWvattl 
2005) or the sole interaction between the dust and the gas 
dTakeuchi & Artvmowiczl200ll) . In the two cases, the authors do 
not consider the influence of the two companion stars which can 
be justified by our study of the kinematic constraints. It is indeed 
possible to obtain a binary orbit reproducing the observed radial 
velocities but which does not let the two M stars generate any 
perturbation to the disk because they pass too far from the disk 
or beca use they have not yet passed at the periapsis. The planet 
used by IWvattI (l2005h is efficient to open the gap and generate 
the outer ring but, as shown in Sect. |5] it is difficult to improve 
this scenario i n order to also explain the inner r ing and the inner 
disk depletion. iTakeuchi & Artvmowiczl(l2001h . with an alterna- 
tive scenario of dust migration due to the gas friction, generate 
a two ring structures but the hypothesis of axysmetricity does 
not leave a lot of possibilities to extend this work and explain 
the tightly-wound spiral structure of the outer ring nor the spiral 
ar m in the extended diff use emission. 

lArdila et al.l ( l2005h consider several physical processes for 
the modelling of HD 141569: gravitation, interaction between 
the gas and the dust, dust generation by collisions. It is there- 
fore, up to now, the most complete theoretical study about HD 
141569. They considered the case of a flyby in order to obtain 
large spiral arms similar to those observed in the extended dif- 
fuse emission. They also added a 5 My planet, on an orbit with a 
semi-major axis of 100 AU with an eccentricity of 0.6, to explain 
the depleted zone below 150 AU. Altought the authors obtain 
good results for the the extended spiral arms and the inner hole, 
they did not succeed to reproduce the two rings nor the gap in be- 
tween. One possible explanation is that, in our work, we obtain 
satisfactory results only if we let the disk evolves with the planet 
during several million years before the flyby. To go further in the 
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description of the HD 141569system, it will be therefore inter- 
esting to use the model developed by these authors but with the 
initial conditions on the binary orbit and on the planet parame- 
ters established by our work. 

7. Summary and conclusion 

The binary orbit is finally found to be almost fixed by the ob- 
servational constraint on a edge-on plane with respect to the ob- 
server. If the binary has had an influence on the disk structure, 
it should have a passing time at the periapsis between 5000 and 
8000 years ago and a distance at periapsis between 600 and 900 
AU. The scenarios with retrograde orbits better reproduce the 
observations but prograde orbits cannot be totally excluded. In 
order to reproduce the observed structures in the debris disk, it 
also appears that the best scenario is a flyby with 1 planet embed- 
ded in the disk. For a 2 My planet, its orbital eccentricity must 
be around 0.2 while for a 8 My planet, it must be below 0.1. In 
the two cases, its apoapsis is about 130 AU. 

This scenario is able to both reproduce the observed radial 
velocities of the stars and the overall structure of the disk with 
the modelling of a flyby by the binary and a embedded planet 
in the disk. However, this scenario is not fully satisfactory be- 
cause it does not reproduce all the structures, and in particular 
the entire gap between the two annulus. Future study can there- 
fore extend this scenario, for example with the modelling of the 
gas and/or collisions, as we have not taken into account these 
efl'ects on the dust dynamics. 
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